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FORM FACTORS OF COMPOSITE SYSTEMS 
BY GENERALIZED WIGNER-ECKART THEOREM FOR 

POINCARE GROUP 

The relativistic approach to electroweak properties of two-particle composite systems 
developed previously is generalized here to the case of nonzero spin. This approach is based 
on the instant form of relativistic Hamiltonian dynamics. A special mathematical technique 
is used for the parametrization of matrix elements of electroweak current operators in terms of 
form factors. The parametrization is a realization of the generalized Wigner-Eckart theorem 
on the Poincare group, form factors are corresponding reduced matrix elements and they have 
the sense of distributions (generalized functions). The electroweak current matrix element 
satisfies the relativistic covariance conditions and in the case of electromagnetic current it 
q , also automatically satisfies the conservation law. 
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1 Introduction 



The development of correct quantitative methods of calculation of composite-particle struc- 
ture is an important line of investigations in particle physics. The modern experimental 
researches of the particle structure are performed with such an accuracy and for such en- 
ergies that it calls for the relativistic theoretical method of interpretation of results. The 
important feature of these methods is the extraction of kinematical parts and invariant parts 
(form factors) from the current matrix elements. A representation of the matrix element in 
terms of form factors is called the parametrization of the current operator matrix element. 
Let us remark that it is just form factors - Lorentz invariant functions - that are extracted 
from scattering data usually. 

In the relativistic theory of the description of composite systems, it is possible to identify 
two main but very different approaches (see e.g. and the pictorial scheme on Fig. 11 in 
the review 

The first is the method of field theory. Based on the principles of quantum field theory - 
quantum chromodynamics (QCD) - it is rightly regarded as the most consistent approach to 
the solution of this problem. However, standard perturbative QCD gives sufficiently reliable 
computational prescriptions only for the description of so-called "hard" processes, which 
are characterized by large momentum transfers, and it does not permit the calculation of 
characteristics determined by "soft" processes. Moreover, there are strong indications that 
perturbative QCD is not valid for the description of the currently existing experimental 
facts in exclusive processes. This applies, in particular, to the description of the elastic form 
factors of such well-studied composite systems as pion, nucleon and deuteron. 

The second method in the relativistic theory of composite systems that is gaining accep- 
tance in recent years is based on the direct realization of the algebra of the Poincare group 
on the set of dynamical observables on the Hilbert state space of the system. This approach 
is called the theory of direct interaction, or relativistic Hamiltonian dynamics (RHD). The 
discussion of this approach can be found e.g. in [21 El- 

The method is based on the quite different starting positions than the first one. It is one 
of the reasons why the results of these approaches are difficult to compare. We shall not 
focus ourselves on the details of RHD but we refer reader to the paper [I] and references 
therein. 

Although the above-mentioned methods of the description of composite systems are dif- 
ferent enough, it is necessary to solve in both cases three highly important problems. 

1) It is necessary to give the relativistic recipe for taking into account the nonperturbative 
effects. Notice that in the case of RHD, in the framework of which we shall operate, these 
effects are contained in the structure of method directly [TJ H] . 

2) It is necessary to construct the matrix element of electroweak current with correct 
transformation properties and to perform the procedure of the relativistic parametrization 
of matrix elements that is to extract the relativistic invariant form factors. 

3) It is necessary to take into account the fact that the form factors of composite sys- 
tems have, generally speaking, the sense of distributions(generalized functions), even in the 
nonrelativistic case. 

In simple cases the parametrization can be obtained with reasonable facility from general 
considerations. However, in more complicated cases, for example, in the case of composite 
systems with arbitrary total angular momentum one needs a special mathematical method. 

The aim of the present paper is the construction of a general parametrization method for 
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the case of composite systems with nonzero values of total spin. 

The method was announced in our brief report [Sj. Now we shall give the consistent 
presentation of our method which makes it possible to apply it to the specific composite 
systems. Our method is similar in some sense to the method of canonical parametrization of 
one-particle matrix elements of field operators It should be particularly emphasized that 
we are dealing now with relativistic Hamiltonian dynamics rather than with field theory as 
inRef. 0. 

Conceptually, our method of parametrization is a realization of the Wigner-Eckart theo- 
rem for the Poincare group and form factors are the reduced matrix elements. This theorem 
is used extensively in the different fields of theoretical and mathematical physics from well 
known theory of angular momentum to the Hopf algebras and quantum superalgebras [Zj. 

In the present paper we use the term - generalized Wigner-Eckart theorem for the 
Poincare group - when form factors of composite systems are considered as generalized 
functions. The necessity of such consideration was shown in Refs. pi E] in the study of form 
factors of composite systems with zero total spin. 

The paper is organized as follows. 

In Sec. II the canonical parametrization of local operator matrix elements between one- 
particle states of arbitrary nonzero spin is described. This parametrization presents the 
extraction of the reduced matrix element on the Poincare group. The electromagnetic current 
matrix element is derived for the spin 1/2 particle. 

In Sec. Ill we show how to construct the electromagnetic current operator for composite 
system of two free particles of spin 1/2. The current matrix element is obtained in the basis 
where the center-of-mass motion is separated. The electromagnetic properties of the system 
are defined by reduced matrix elements, or the so called free two-particle form factors, these 
form factors being generalized functions. This means that, for example, the static properties 
of the system are given by the weak limits as Q 2 — > (Q 2 = —q 2 , q is the momentum 
transfer). A special attention is paid to the case of total spin one and zero total orbital 
momentum. In this case the electromagnetic properties are defined by four free two-particle 
form factors - charge, quadrupole, magnetic and the magnetic quadrupole form factor of the 
second kind (see the review [H] and the references therein). 

In Sec. IV the procedure of the construction of the electromagnetic current matrix ele- 
ment is developed for the case of two interacting particles. Each step of the calculation in 
the developed method remains Lorentz covariant and conserves the current. The compos- 
ite system form factors are derived as some integral representations in terms of the wave 
functions obtained in the frame of RHD. 

In Sec.V the results are summarized and the conclusion is given. 

2 Parametrization of one— particle operator matrix el- 
ements 

Let us describe now a general method of canonical parameterization of local operator matrix 
elements. In the present paper we are dealing with the matrix elements diagonal in total 
angular momentum only. 

The main idea of the parametrization can be formulated as follows. Using the variables 
entering the state vectors which define the matrix elements one has to construct two types 
of objects. 
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1. A set of linearly independent matrices which are Lorentz scalars (scalars or pseu- 
doscalars). This set describes transition matrix elements non-diagonal in spin projections 
in the initial and finite states, as well as the properties defined by the discrete space-time 
transformations. 

2. A set of linearly independent objects with the same tensor dimension as the opera- 
tor under consideration (for example, four-vector, or four-tensor of some rank). This set 
describes the matrix element behavior under the action of Lorentz group transformations. 

The operator matrix element is written as a sum of all possible objects of the first type 
multiplied by all possible objects of the second type. The coefficients in this representation 
as a sum are just the reduced matrix elements - form factors. 

The obtained representation is then modified with the use of additional conditions for 
the operator, such as the conservation laws, for example. In order to satisfy these additional 
conditions in some cases some of the coefficients (form factors) occur to be zero. 

To demonstrate this let us consider the parameterization of the matrix elements taken 
between the states of a free particle of mass M in different simple cases. Let us normalize 
the state vectors as follows: 

(p,m\p' ,m') = 2p 5(p- p') S mm > , (1) 

here p , p' are three-momenta of particle,po — \IM 2 + p 2 , m , m' are spin projections. 

Let us consider first the parameterization of the matrix element of scalar operator A(x) 
taken between the states of free particle with zero spin. Because of translation invariance it 
is sufficient to consider A(0). Using the variables of state vectors in the initial and final states 
one can construct two linearly independent scalars: p 2 = p' 2 = M 2 , {p — p') 2 = q 2 = —Q 2 . 
Only the second one - the momentum transfer square - is nontrivial. Ignoring the trivial 
dependence on M we can write: 

(p, M\A(0)\p',M) = f(Q 2 ). (2) 

Let us consider now the matrix element of a scalar operator between the states of a 
particle with spin j: 

(p,M,j,m\A{0)\p',M,j,m'). (3) 

The tensor dimension of the operator is not changed, nevertheless, the matrix element is 
now a matrix in spin projections in the initial and finite states. 

Let us construct the set of Lorentz scalars - linearly independent matrices in spin pro- 
jections to be used for the construction of the representation of Eq. Q. Let us use the 
covariant spin operator T M (p). In the rest frame this operator coincides with the particle 
spin operator: 

r°(0) = 0, f(0)=Mj, [j i J h ]=ie m j l . (4) 

The covariant spin operator can be defined with the use of the Pauli-Lubanski vector u> M [Oj. 
In terms of matrix elements we have: 

^ram'ip) = (P ™KIP , ™'> = < , H U~ X (A>^(A p ) |0 , rri) = 

= (AX(0,rn\w»\0,m') = (A p )£ 1^(0) . (5) 

Here A p is the boost for the transformation from the rest frame of the particle to the lab- 
oratory frame, U(A P ) is the corresponding representation operator. The matrix A p can be 
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written explicitly using the matrix for Lorentz transformation of the vector p' into the vector 

A.^ + ^-^M. (6) 
" v M* P Pv M 2 +p x p' x K } 

In our case A p is given by Eq. (JHJ) with p' = (M , , , 0). 
Using the explicit form of A p we obtain: 

T (p) = (pj), f{p) = M 3 + ^\- r 2 = -M 2 j(j + l). (7) 

PO + M 

As is known JH], spin transforms under the action of Lorentz group following the so 
called little group which is isomorphic to the rotation group, that is the corresponding 
transformations are realized by the matrices of three dimensional rotations. 

So, the 4-spin operator is transformed under Lorentz transformations p^ = Af\, p' v in the 
following way: 

r"(p) = K d\p, p') rv) D\ P \ P ) . (8) 

Using Eq. (JHJ) one can show directly that the matrix elements of the operators D^{p, p')T^{p') 
and T fJ, (p)D : ' (p, p') transformed as 4-pseudovectors, the matrix elements of the operators 
Di(p, p')p fl T tJ '(p') and p' r M (p)I? J '(p, p') - as 4-pseudoscalars. 

We shall construct the set of linearly independent Lorentz-scalar matrices using the vec- 
tors p ,fl and the pseudovector Z? 3 '(p, p')Y^{p'). Note, that the pseudovector r M (p)D J (p, p') 
does not enter the decomposition, being linearly dependent. One can show this fact using 
the relation (JHJ) and the explicit form Eq. (JHJ) of A^, transforming p' into p. It is easy to 
obtain: 

p^ + p ,fl 
M 2 + p^p 1 ' 

As p'^T^lp') = 0, the set in question of linearly independent matrices in spin projections of 
the initial and the final states giving the set of independent Lorentz scalars is presented by 
2j + 1 quantities: 

D j (p, p') {p^{p')) n , n = 0,l,...,2j. (10) 

The number of linearly independent scalars in Eq. (fTU|) is limited by the fact that the product 
of more than 2j elements r M (p') reduces as it is linearly dependent. If n is even, then the 
obtained quantities are scalars, if n is odd then they are pseudoscalars. The current matrix 
element (JHJ) is represented by the linear combination of these linearly independent Lorenrz 
scalars. The coefficients in this combination f n {Q 2 ) are just form factors. These form factors 
are invariant under rotations and so they do not depend on spin projections. So, they depend 
upon only one scalar combination of variables - the momentum-transfer squared. 

For self-adjoint operator A(0) a minor modification of the set ()10)1 is necessary: in the 
scalar product in Eq. (JTH)) the factor i appears. So, now the current matrix element (JHJ) can 
be written in the form: 

(p,M,j,m\A(0)\p r , M,j, m') = 
= E E (m\D^p,p')\m'')(m''\{ l p^(p')r\m')f n (Q 2 ) (11) 

ra=0 m"=—j 

The obtained operator is self-adjoint. This can be easily shown using the following relation: 

p f r»(p)Di(p, p') = -D*(p, p')p^{p') (12) 



r»z^, p>) = D*(p, p>) {rv) [p„rV)] • (9) 



5 



In the case of a scalar operator the values n will be even and for a pseudoscalar they will be 
odd. 

Let us consider now the 4- vector operator ^(0). To parametrize the matrix element one 
needs a set of quantities of the appropriate tensor dimension. Using the variables entering 
the particle state vectors one can construct one pseudovector r M (p') and three independent 
vectors: 

K, = (p-p'), = q,, #; = (p + 2>V R, = e, uXp p"p lX V{p') . (13) 

Here t pv \ p is a completely anti-symmetric pseudo-tensor in four dimensional space-time with 
^0123 — — 1- 

The operator matrix element contains the matrix elements of the listed quantities mul- 
tiplied by D^(p, p') from the left. Each of such products is to be multiplied by the sum of 
linearly independent scalars constructed while obtaining the parameterization !)!!)) : 

(p, M, j, m |j M (0) | p', M,j, m') = 

= J2 ( m\D>(p, p') \m") ( m"\ F 1 K'^ + F 2 T^p') + F 3 R P + F 4 K^m') , (14) 

to" 

where 

*i = £ fin{Q 2 )(ip^{p')Y ■ (15) 

n=0 

Let us impose some additional conditions on the operator. 

1. Let us require the operator to be self-adjoint. 

One can check this condition by use of Eq. (|12jh In the right-hand side (r.h.s.) of 
Eq. we need to modify slightly the vector multiplied by F 2 . The new vector is a linear 
combination of the 4-vectors entering Eq. ()14jl and has the following form: 

r<V) - rv) - §^ (pjrtf)) ■ (16) 

The terms containing F 2 and F3 are modified in the following way: 

FiA» -> ^(F t A» + A»F l ) = {F t A»} + , 2 = 2,3. (17) 

Here N 1 are the vectors entering Eq. ((H)) and changed following Eq. (jlfij) . The terms 
containing F 3 and _F 4 have to be multiplied by i. 

2. It is useful to modify the parametrization so as to make the vectors - multipliers of 
Fi - orthogonal to each other. 

This results in the following form of the vector multiplied by F 2 ([IS)) : 

(K'* 

rV) -[-^2+^) (p^CpO) • (is) 

3. Let us impose the condition of parity conservation. 

This condition is satisfied if all the terms in the sum (JT5j) contain an even number of 
pseudovector factors T M . So, the sums in F\ and F4 are over even n : 2j > n > 0, in F3 
over even n : 2j — 1 > n > 0, and in F 2 over odd n : 2j — 1 > n > 0. 

4. Let us impose the conservation condition j^K^ = j p q^ = 0. 
It is easy to see that this condition is satisfied only if F 4 = 0. 
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So, finally we have the following parametrization of the matrix element: 
(p, M, j, m 1^(0)1 p', M, j, rri) =Y,i™\D j (P, p')\m")x 



x ( m"\ + {F 2 -(p^(p')) (§2 + ^ 



+ i{F 3 RA,\m'), (19) 



This construction can be used, for example, to obtain the electromagnetic current matrix 
element in the case of particle with spin 1/2. 

Now let us list the conditions for the electromagnetic current operator to be fulfilled in 
relativistic case (see Ref. [3] and the references therein). 

1. Lorentz covariance: 

U-\A)f(x)U(A) = KfiA-'x) . (20) 

Here A is a Lorentz-transformation matrix, U (A) is an operator of the unitary representation 
of the Lorentz group. 

2. Invariance under translation: 

U-\a)nx)U(a)=f(x-a). (21) 
Here U (a) is an operator of unitary representation of the translation group. 

3. Current conservation law: 

[P u f(0)} = 0. (22) 
In terms of matrix elements ( j^ 1 ) the conservation law can be written in the following form: 

ff„<?(0)> = 0. (23) 

Here is four-vector of the momentum transfer. 

4. Current-operator transformations under space-time reflections: 

U P ^j\x\x)3{x\x)jUp l = (jV,-f),-J(x°,-f)) , 

U R y{x)Un X =f{-x). (24) 
In Eq. (I24J) Up is the unitary operator for the representation of space reflections and Ur is 



the anti-unitary operator of the representation of space-time reflections R = PT. 

Our parametrization satisfies all the conditions listed above. The parametrization fTH|) 
for the electromagnetic current in the case of particle with spin 1/2 has the following form: 

(p,M, -, m|j M (0)|p', M, \,m') = 
= Y,(™\D^(p, p')\m")(m"\ f 10 (Q 2 ) + i f 30 (Q 2 ) R,\m>) , (25) 

m" 

The form factors fw(Q 2 ) = fi(Q 2 ) and fw(Q 2 ) = f2(Q 2 ) are the electric and the magnetic 
form factors of the particle, respectively. In analogous way (although a little more cumber- 
some) one can obtain the matrix elements of the operators of higher tensor dimension. 
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3 Parametrization of matrix elements of two— particle 
electromagnetic current operator 

To describe the properties of the system of interacting constituents in our approach it is 
necessary to have the reduced matrix elements on Poincare group (the form factors) which 
describe the properties of the composite system of free constituents. In this section we 
generalize the method of parametrization of the previous section to the case of such free 
systems. 

Let us consider a system of two free particles with spins 1/2 and let us parametrize the 
matrix element describing the transitions in this system. Let us construct, for example, the 
electromagnetic current operator matrix element. The matrix element can be taken between 
the following two-particle state vectors: 

| pi , mi, p 2 , m 2 ) = | pi , mi) <g> |p 2 , m 2 ) . (26) 

Here pi , p 2 are three-momenta of particles, ra\ , m 2 are spin projections on the z axis. The 
one-particle state vectors are normalized by Eq. (P). 

As well as the basis (|2fij) one can choose another set of two-particle state vectors where 
the motion of the two-particle center of mass is separated: 

| P, y/s, J, I, S, mj) , 

(P, \/s, J, I, S, mj\P', vV, J', I', S', m j') = 

= N CG 6®(P - P')5{^Ts - \fa)6jj,6 U r6ss<6 mjmjl , (27) 



(2P ) 2 . y/\(s,MZ,M2: 

iv CG — — — ■, & 



here P M = (pi + p 2 )^, Pf L = s, yfs is the invariant mass of the two-particle system, / is the 
orbital angular momentum in the center-of-mass frame (cm.), S is the total spin in the 
cm., J is the total angular momentum with the projection mj; M\ , M 2 are the constituent 
masses, and A(a, b, c) = a 2 + b 2 + c 2 — 2(ab + be + ac). 

The basis (|2*Tj) is connected with the basis (j2El) through the Clebsh-Gordan decomposition 
for the Poincare group. Here we write the decomposition in a little more general form than 
in Ref. gj: 

| P, y/s, J, I, S, mj) = 
I ^ ^ IPi > m i; P2,m 2 )(p 1 ,m 1 ; p 2} m 2 \P } y/s, J, I, S, mj) , (28) 



mi ma 

where 



[p 1 ,m 1 ;p 2 ,m 2 \P, y/7s, J, I, S, mj ) = V2~s[\(s, M 2 , M 2 )]" 1 / 2 2P 5(P - pi - p 2 ) x 
x J2( m i\ P> 1/2 (Pi , P) \mi)(m 2 \ D 1/2 (p 2 , P) \m 2 ) x 
x(- ^mim 2 \Sm s )Yi mi ({} ,ip) (Slm s mi \Jmj) . 
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Here $ , <p are the spherical angles of the vector p = (pi — P2V2 in the cm., F; mi ($, <p) 
are spherical harmonics, (1/2 1/2 mi m 2 \Sms ) , (S lm s mi\Jmj) are Clebsh-Gordan co- 
efficients for the group 577(2), U J are the known rotation matrices to be used for correct 
relativistic invariant spin addition, the summation is over mi, m 2 , mi,ms ■ 

The decomposition in spherical harmonics and angular momenta summation in Eq. (|28jl 
are performed in the cm. and the result is shifted to an arbitrary frame by use of D- 
functions. 

The electromagnetic current matrix element for the system of two free particles taken in 
the basis (|2fij) can be written as a sum of the one-particle current matrix elements: 

(pi, mi;p 2 , m 2 \ jj?\0) \p[, mi;p 2 , m 2 ) = 

= (P2,m 2 |p 2 ,m 2 )(pi,m 1 | < 7'i M (0)|p / 1 ,m / 1 ) + (1^2). (29) 

Each of the one-particle current matrix elements in Eq. (|2*9*j) can be written in terms of form 
factors as in Section II. In the case of the particles with spin 1/2 we make use of Eq.(|25J). 
So, in this case the electromagnetic properties of the system are defined by the form factors 

Now let us construct the electromagnetic current matrix element for the system of two 
free particles in the basis ()27|) following the previous Section. Let us consider first a simple 
case J = J' = S = S' = I = I' = 0. (We omit these variables in the state vectors.) 
This set of quantum numbers appears, for example, in the case of pion Now there is 
no pseudovector T^, but along with the scalar (P — P') 2 = — Q 2 two additional nontrivial 
scalars do appear s' = P' 2 and s = P 2 - the invariant mass squares for the free two-particle 
system in the initial and in the final states. So, the form factors entering the parametrization 
are functions of the variables Q 2 , s , s'. The current matrix element is presented by the 
linear combination of the four- vectors P.. and P'\ 

P [A 

(P, V~s I jf (0) I P', = (P„ + Pi) gi(s, Q 2 , s') + (P, - Pi) g 2 (s, Q 2 , s') . (30) 

Making use of the conservation condition (|2Hj) 

4°)(0)(P '- P'Y = , (31) 

we can write the parametrization in the form: 

(P, yT s I ^)( ) I P> 7 y/P) = A ,(s, Q 2 , s') g (s, Q 2 , s') . (32) 

Here go(s, Q 2 , s') is the reduced matrix element. We will refer to this invariant as to the free 
two-particle form factor. The vector A^(s, Q 2 , s') is defined by the current transformation 
properties (the Lorentz covariance and the conservation law): 

A,(s, Q 2 , s') = (1/Q 2 )[(s -s' + Q 2 )P, + {s'-s + Q 2 )P'^} . (33) 

The free two-particle form factor can be expressed in terms of the one-particle form 
factors (|25j) . To do this one has to perform in Eq. (}32|) the Clebsh-Gordan decomposition of 
the irreducible representation ()27j) into the direct product of two irreducible representations 
(J2*HJ) . (|2*K|) and to take into account Eqs. (j25)l . (|29)1. As the form factors are invariants one 
can perform the integration in Eq. (}32|) in the coordinate frame with P' = , P = (0, 0, P). 
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The explicit form of go(s, Q 2 , s') in the case of two particles with spin 1/2 and mass M can 
be found in Ref. jlj. 

Let us perform the analogous parameterization for the set of quantum numbers in the 
basis (J2*7)) in the case J , J' ^ 0. In the following, we will take J = J' . 

The Lorentz covariant properties of the matrix element are defined (in analogy with 
Eq. (JUDO by three 4-vectors and one pseudovector: 



K'=(P + P% , K, = {P- P% , R, = 6^ P P V P' X TP{P') , T M (P') 



(34) 



The pseudovector r M (P) does not enter the parametrization because it can be expressed 
through r^(P') by the equation analogous to Eq. Q: 



-= r(p)„ d'(p, p') = d'(p, p') -=r„(P 



ss' P V P' v + Vss 

The set of linearly independent matrices to be used for the decomposition of the current 
matrix element is obtained from the vectors P^ and T /i (P / ) following Eq. (fTU|) : 



D J (P, P') ( P M r^(P')) n , n = 0,1,..., 2 J 



(35) 



Using the conditions of self-adjointness, current conservation ()31|) . parity conservation 
and orthogonality of the parametrization vectors to one another (as in Eq. (fT9"jO we obtain 
(see Eq. (|T7)l. too): 

(P, v^, J, Z, 5, |P', vV, J, Z', 5', mS) = 



^(m J |^(P,P')|m'i)(m';| £ {*f S5 ' 4( a » <?V)}_ 

m'J 

1 



(36) 



A 1 



~ Q 2 



s-s' + Q 2 )P, + {s'-s + Q 2 )P' 



2v^i 



A( s ,-g 2 , S ') 



pp' 

P\T\P') 



A 



Pn 



The quantities Ff ss ' in Eq. (jHEJ) are defined by the relations analogous to Eq. ()15j) : 



2./ 



^/Z'SS" 



E^ 5S '(^Q 2 ^')(^(Por. 



(37) 



(3? 



n=0 



The sum in Eq. (J38J) is taken using the parity conservation condition as in Eq. (|T9"j) . 

Let us remark that the reduced matrix elements - the invariant form factors — now (in 
contrast with the form factor ()32)l ) depend on the additional invariant quantities I, I', S, S' 
that are invariant degeneration parameters in the basis ()27j) . 

The self-adjointness condition (}3l)j) is fulfilled at following constraints on form factors 

ll'SS'. 



f, 



f 



ll'SS' 



s,Q 2 ,s')=fi;r s (s>,Q 2 ,s). 



(39) 
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Here the star means the complex conjugation. 

Let us consider especially the case M 1 = M 2 = M , J = J' = S = S' =1 , I = I' =0. 
The reduced matrix elements for this case can be used to calculate p - meson properties 
neglecting the D - state contribution. In this case the functions Ff ss in Eq. (f3*K|) have the 
following form (compare to Eqs. (fTHjl . (fTHjl ): 

Fx = Ma, Q 2 , s') + f 12 (s, Q 2 , s')(iP v T\P')) 2 , 

F 2 = f 21 {s,Q 2 ,s'){iP v T v {P')), F 3 = f 30 (s,Q 2 ,s'). (40) 

In equations (|4U|) the fixed variables I, I', S, S' are omitted. Time reflection invariance imposes 
the additional conditions: 

f* n = f m , i = / 2 * 1 = -/ 2 i. (41) 

With Eqs. (JHJ) and (jUJ) we obtain 

f in (s,Q 2 ,s') = fUs',Q 2 ,s), z = l,3; f2i(s,Q 2 ,s') = -f 21 (s',Q 2 ,s). (42) 



The relation (}42j) demonstrates that the form factor / 21 appears in the parametrization as 
a consequence of the fact that the invariant masses in the initial and the final states are 
different. As will be seen below, this form factor gives no contribution to elastic processes, 
for example, to the electron scattering by the composite particle, however, it does contribute 
to radiative transitions. 

Let us rewrite Eq. (j4T)j) in terms of standard form factors instead of fi n : 

F, = goc{s, Q 2 , s') + g 0Q (s, Q 2 , s') {(zP^(P')) 2 -(1/3) Sp(zP^(P')) 2 } gp(p ^ (p/))2 , 

F 2 = g 0MQ (s, Q 2 , s') (i P V T V (P')) , F 3 = g 0M (s, Q 2 , s') . (43) 

The scalar factor in the term with g Q is chosen in such a way as to make it possible to 
interpret g$Q as a quadrupole form factor of the system of two free particles. In other terms 
goc is the charge form factor, goMQ = /21 is the magnetic quadrupole form factor of the 
second kind, its classical analog being the so called toroidal magnetic moment |H] , goM = /30 
is the magnetic form factor of the free two- particle system. 

In the same way as in the case of Eq. (}3*2*|) the free two-particle form factors in Eq. 
can be written in terms of the constituent form factors (see Ref. |3] for details). The 
corresponding equations are rather complicated, so we do not present it here. 

The free two-particle form factors in Eqs. (|32)1 . (j36j) . (|4()|) are to be considered in the 
sense of distributions. For example, go(s , Q 2 , s') has to be interpreted as a Lorentz invariant 
regular generalized function on the space of test functions <S(R 2 ) [TU] . 

Let us define the functional giving the regular generalized function as 

( g (s, Q 2 , s') , <p(s, s')) = J dp(s, s') g (s, Q 2 , s') ip{s, s') . (44) 

Here 

dp{s,s') = lQd(s-AM 2 )d(s'-AM 2 )<f^'dp{s)dp{s') , dp(s) = -k d^s . (45) 
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The quantity Q 2 is a parameter of the generalized function, Mi = M 2 = M. The function 
9{x) is the step function. 

if(s , s') is a function from the space of test functions. So, for example, the limit as 
Q 2 — > (the static limit) in g (s , Q 2 , s') (jH2j) exists only in the weak sense as the limit of 
the functional: 

lim (g ,<p) = ((e q + eq)6(n(s') - fi(s)),(p) . (46) 

e g and e 5 are the constituent charges, 5 is the Dirac delta-function. 

It is just the limit in the sense (|46j) that gives the electric charge of the free two-particle 
system. The ordinary point-wise limit of the form factor g as Q 2 — *>0 is zero. The equations 
analogous to Eq. ()46|) are valid for the static limits of the free two-particle form factors in 
Eqs. (JliJ), too. 

Let us note that the conditions imposed on the free two-particle form factors that follow 
from the conditions for the electromagnetic current operator . and (|42|) have to be 
considered in the weak sense, too. 

4 Parameterization of the current operator matrix el- 
ements for systems of two interacting particles 

In this Section we generalize the parameterization method of the previous sections to the 
case of composite system with the structure defined by the interaction of its constituents. 

Let us consider the operator j M (0) that describes a transition between two states of a 
composite two- constituent system, j M (0) being diagonal in the total angular momentum. 
Let us neglect temporarily the additional conditions of self-adjointness, parity conservation 
etc. in the same way as when constructing the matrix elements (|T4"j) in Section II. The 
Wigner-Eckart decomposition of the matrix element has the form (|14jl . |T5|) . To emphasize 
the fact that the particle is composite, let us rewrite Eqs. Ifl4*jl . (fTKJ) using new notations: 

(Pc m Jc |j M (0) | p c ',rn'j c ) = 

= (m Jc \D J c( Pc , p' c ) [F?Kl + F°TM +F^R, + F^K,] \w! Jc ) , (47) 

here 

2J C 

F° = Tftn{Q 2 )(iVc^{p' c )Y. (48) 

n=0 

In Eqs. p7j) . and (|4*Hjl {p c —p'c) 2 = — Q 2 , P%= Pcu 2 = M 2 , M c is the mass of the composite 
particle. In contrast to Eq. (|T4*j) spins and masses are omitted in the state vector variables. 

In the frame of RHD the form factors of composite systems f£ n are to be expressed in 
terms of RHD wave functions 

In RHD a state of two particle interacting system is described by a vector in the direct 
product of two one-particle Hilbert spaces (see, e.g., Ref. |lj). So, the matrix element in 
RHD can be decomposed in the basis (|27jl : 

{Pc,mj c \j^)\p c ' ,m' Jc ) = 

f dPdP' 

/ t^— dy/s dy/s* (p c ,m Jc \P ,y/s ,J ,1 ,S ,mj)x 

J IV CG ^CG 
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x{P,^,J ,l,S ,mj\j^O)\P' ,f ,1' ,S' , mj/ )x 

x (P',VP, J' ,1' ,S' , mji \p c ' , m' Jc ) . (49) 

Here the sum is over variables J,J',l,l',S,S',mj,mji, and (P ' , yfs 1 J' \V ', 5" ,m'j\p c ' ,m' Jc ) is 
the wave function in the sense of the instant form of RHD. 

Let us write the wave function in the form slightly more general than in Ref. [3]: 

(P , \/s ,J ,1 ,S ,mj\p c ,m Jc ) = N C 5(P - Pc)Sj c j5 mjcmj ^§( k ) ■ ( 50 ) 




The RHD wave function of constituents with fixed total angular momentum and total spin 
is defined by the form: 

rf§(k(s)) = y/Vt(l-rj*/s*)u ls (k) k , (51) 
and is normalized by the condition: 

£ / uUk)k 2 dk = l. (52) 
is 

Here r\ = Mf — Mf ,uis(k) is a model phenomenological wave function. 

To calculate the form factors in Eqs. (|47j) . and (|48jl let us write the Wigner-Eckart 
decomposition on the Poincare group for the current matrix element in the r.h.s. of Eq. ()49|) . 
However, now there are some difficulties. 

The point is that in the previous sections we were dealing with the parametrization of 
local operator matrix elements in the case when the transformations of the state vectors and 
of the operators were defined by one and the same representation of the quantum mechanical 
Poincare group. While describing the composite systems in RHD a different situation can 
arise when the state vectors and the operator under consideration are transformed following 
different representations of this group. 

It is just such a situation that takes place in the case of the matrix element in the r.h.s. of 
Eq. ()49j) . The operator in the matrix element describes the system of two interacting particles 
and transforms following the representation with Lorentz boosts generators depending on 
the interaction [3]. The state vectors physically describe the system of two free particles 
and present the basis of a representation with interaction-independent generators. So, the 
Wigner-Eckart decomposition in the form used in Section II can not be applied directly to 
the matrix element in the integrand in the r.h.s. of Eq. (}49|) . This is caused by the fact that it 
is impossible to construct 4-vectors describing the matrix element transformation properties 
under the action of Lorentz boosts from the variables entering the state vectors (contrary 
to the case of, e.g., Eq. ()47|) ) . In fact, the possibility of matrix element representation in 
the form (}4Tj) is based on the following fact. Let us act by Lorentz transformation on the 
operator: 

f>- 1 (A)^(0)f>(A)=?(0). (53) 
We obtain the following chain of equalities: 

(p c ,m Jc \f(0)\p c ',m' Jc ) = (p c ,m Jc \U-\A)f(0)U{A)\p c ', m ' Jc } = 
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= £ (m Jc \[D J '{R K )]- x \m Jc )(kp^ (54) 

m Jc ,m' Jc 

Here D Jc (R A ) is rotation matrix realizing the angular momentum transformation under 
the action of Lorentz transformations. The equalities (p>fj) show that the transformation 
properties of the current as a 4-vector (|20|) can be described using the 4-vectors of the 
initial and the final states. This means that the procedure of the canonical parameterization 
can be extended to the matrix element for the realization of the Wigner-Eckart theorem on 
the Poincare group. 

In the case of the current matrix element in the r.h.s. of Eq. ()49|) the relations (p>4*|) are 
not valid and direct application of the Wigner-Eckart theorem is impossible. 

However, it can be shown that for the matrix element in Eq. (J4*§|) considered as a gener- 
alized function (that is considered as an object having sense only under integrals and sums 
in Eq. (JHUO, the equality (jHljl is valid in the weak sense. 

Let us consider the matrix element in question as a regular Lorentz covariant generalized 
function (see, e.g., Ref. [TU] ) . Using Eq. (joUj) let us rewrite Eq. (|4^j) in the following form: 

(pc,mj e \j IM {Q)\p e , ,m' JC )= E / AT c ' dyfidVtf<pfS(s)<Pi. c s>(s')x 

l,l',S,S' J ^CG^CG 

x (p c ,y/s, J c ,l, S ,m Jc |j M (0)|p c ', v 7 ^ 7 , J c ,/', 5", m' Jc ) . (55) 

Here it is taken into account that the current operator j M (0) is diagonal in total angular 
momentum of the composite system. 

Let us make use of the fact that the set of the states ([27)1 is complete in RHD: 

/dP -> -> 

— — d<Js\P ,>/s,J ,l,S ,mj)(P ,yfs,J ,l,S ,mj\ . (56) 
^CG 

Here the sum is over all the discrete variables of the basis (|27jl. 

Under the integral the matrix element of the transformed current satisfies the following 
chain of equalities (flHOJ) and (JoTlj) are taken into account): 

r MM' 

E / N C N : dv~sdv?^(s)^us')x 
j i\ CG i\ CG 

x(p c ,^,J C} l,S, m Jc \U-\k)3^)U(K)\p c ' , v 7 ?, J c ,l',S', m' Jc ) = 
= (p c ,mj c \U-\A) I j^O) I U(A)\p c ' ,m' Jc ) = 
= J2 (rn Jc \[D J c{R k )}- l \m Jc )(Ap c ,m Jc \i f\0)I \Ap c ' ,m' Jc ) (m' Jc \D J ^R A )\m' Jc ) = 

rhjc,m' Jc 

r MM' 

= N C N C , dVZd,/s<p%(s)<plMs') E (rn Jc \[D J *{R x )]^\rh Jc )x 

CG CG rh.j c ,fh' Jc 

x(Ap c ,^,J c ,l,S, rhj c \j^0) \Ap c ',V?,J c ,l',S', m' Jc ) ( m' Jc \D J ^R A ) | m' Jc ) . (57) 

It is easy to see that under the integral the current matrix element satisfies the equalities 
analogous to Eq. (jHljl . so now it is possible to use the parameterization method of the 
previous sections under the integral, that is to use the Wigner-Eckart theorem in the weak 
sense. 
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The next step is a parameterization of the matrix element in the r.h.s. of Eq. (|55|). The 
r.h.s. can be written as a functional on the space of test functions of the form (see Eq. (|511 
too)): 

^ l ' ss '(s,s')=u ls (k(s))u l , s ,(k(s')). (58) 

Eq. (J55J) can be rewritten as a functional in R 2 with variables (s, s') (see Eqs. (jHJ)-(Jin|), 

too): 

r N N' 

(Pc,mj c \MO)\Pc',rn' Jc )= E / dn(s,s>) " c ° ^' ss (s,s')x 



x (p c ,y/s ,J c ,l ,S , m Jc \j^(0) \p c ' ,Vs' ,J c ,l' ,S' , m' Jc ) . (59) 

The measure in the integral (|5T?|l is chosen with the account of the relativistic density of 
states, subject to the normalization (Jo"Tj) . (p)2"j) (see Eq. (j^BI too)): 

dfi(s, s 1 ) = 16 6(s - (M 1 + M 2 ) 2 ) 9{s' - (Mi + M 2 ) 2 ) 

x vV^l - v 2 /s 2 )V^(l - V 2 /s' 2 )dfi{s) dfi{s') . (60) 



dfi(s) is given by Eq. (jlBl . 

Notice that the sums over discrete invariant variables can be transformed into integrals 
by introducing the adequate delta-functions. Then the obtained expression is a functional 
in R 6 . 

The functional in the r.h.s. of Eq. ()59|) defines a Lorentz covariant generalized func- 
tion, generated by the current operator matrix element. The integral in Eq. (JoTJj) converges 
uniformly due to the definition (|58j) . 

Taking into account Eq. ()57j) we decompose the matrix element in the r.h.s. of Eq. (J5l?j) 
into the set of linearly independent scalars entering the r.h.s. of Eq. ()47j1 : 

Nc — (p c , y/s, J c , I, S, m Jc \j^0)\p c ' , V?, J c , I', S', m' Jc ) = 



N CG N' CG 

2.1, 

E 

n=0 



m Jc \D J <(p c , p' c ) ^P^(Pc)rK SS '^Q 2 ^')\ m 'jc)^ (61) 



Here A^ (s, Q , s') is a Lorentz covariant generalized function. 

Taking into account Eq. (pTTj) and comparing the r.h.s. of Eq. with Eq. (joTJj) we 
obtain: 

M'.S.S" 

= ( m Jc \ [ f c Xn Kl + f 2n TM +f*n ^ + ftn K»] \m' Jc ) . (62) 

It is easy to see that all the form factors in the r.h.s. of Eq. ()62|) are nonzero if the gener- 
alized function A contains parts that are diagonal (A\) and non-diagonal (A2) in mj c , m' Jc : 

K SS \s, Q 2 , s>) = <f/(s, Q 2 , s>) + Alf(s, Q 2 , s') . (63) 

For the diagonal part we have from Eq. (|fi2"|) : 

E / d^s,s')i; ll ' ss \s,s')(mjMinf'(s,Q 2 ,s')\mj c ) = 

Tin mi J 



1,1' S,S' ' 
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mjc 



fLMK + f c MK» \m Je ). (64) 



The notation ff n [ip] in the r.h.s. emphasizes the fact that form factors of composite systems 
are functionals on the wave functions of the intrinsic motion and so, due to Eq. (|58jl . on the 
test functions. 

Let the equality (JMj) be valid for any test function ip ss (s, s f ). When the test functions 
(the intrinsic motion wave functions) are changed the vectors in the r.h.s. are not changed 
because according to the essence of the parametrization (}4Tj) they do not depend on the 
model for the particle intrinsic structure. So, when the test functions are varied the vector 
of the r.h.s. of Eq. (jMj) remains in the hyperplace defined by the vectors , K'^. 

When test functions are varied arbitrarily the vector in the left-hand side (l.h.s.) of 
Eq. (jMj) can take, in general, an arbitrary direction. So, the requirement of the validity of 
Eq. f!64|) in the whole space of our test functions is that the l.h.s. generalized function should 
have the form: 

Af n f(s, Q\ s>) = Kl Gl ss '(s, Q\ s>) + K, Gl ss '(s, Q\ s>) . (65) 

Here G% ss \s, Q 2 , s') , % = 1,4 are Lorentz invariant generalized functions. 

Substituting Eq. (JfiHj) in Eq. (jo^j) and taking into account Eqs. (JoTj) and (jBHj) we obtain 
the following integral representations: 

fL(Q 2 )= E / dy/5dV?tfs{ a )Gl sa (a.QVte&OO • (66) 

; v g >S r J M!+M 2 

In the case of matrix element in Eq. (|6*2|) non-diagonal in m,j c , m' Jc we can proceed in an 
analogous way and obtain an analogous integral representations for ff n {Q 2 ) , i — 2, 3. 

So, the matrix element in the r.h.s. of Eq. considered as Lorentz covariant generalized 
function can be written as the following decomposition of the type of Wigner-Eckart decom- 
position: 

(jp e ,y/s,J c ,l,S, m Jc \ j M (0) \p e ' , x/s 7 , J c , V , S' , m' Jc ) = 



NcgN'cg 
N C N> 



m Jc \D J c( Pc , p' c ) T x K + T 2 TM +F 3 + ^ K„] \m' Jc ) . (67) 



2J C 

^ = E Gl SS '(s, Q 2 , s'){ip c ^{p' c )Y , (68) 

n=0 

with the constraint (|6*Bj) . 

In Eqs. f!67j) . and (|6*H|l the form factors contain all the information about the physics 
of the transition described by the operator j M (0). They are connected with the composite 
particle form factors (}4Tj) . and PH|) through Eq. ()66|) . In particular, physical approximations 
are formulated in our approach in terms of form factors Gf n ss (s, Q 2 , s') (see Ref. |3] for 
details). The matrix element transformation properties are given by the 4-vectors in the 
r.h.s. of Eq. (JoTj). 

It is worth to emphasize that it is necessary to consider the composite system form factors 
as the functionals generated by the Lorentz invariant generalized functions G l l n ss (s, Q 2 , s'). 

After determination of the vectors A l ^ s {s, t, s') in Eq. (}6*2|) let us impose the additional 
conditions on the matrix elements in Eqs. (}4Tj) . and (j62j) in the same way as we did in Sec. II 
in Eqs. |T6|) - (fT9|) . The r.h.s. of equalities (jlTjl and (jBTj) contain the same 4-vectors and the 
same sets of Lorentz scalars PH|) and (jHEJ), so, to take into account the additional conditions 
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it is necessary to redefine these 4-vectors according to Eqs. (fTfi |) -(fT§ j) . For example, the 
conservation law gives JF 4 =0. It is easy to see that for the redefined form factors the equality 
(I66|) remains valid. 

Let us write the parameterization (|67|) . (fHKJ) for the particular case of composite particle 
electromagnetic current with quantum numbers J — J' — S — S' — 1, which is realized, for 
example, in the case of deuteron. Separating the quadrupole form factor in analogy with 
Eq. and using Eqs. (|S7j). and (J55|) we obtain the following form: 



(p c ,Vs,J c ,l,S, m Jc \ j M (0) \p c ' ,Vs' ,J c ,l' ,S' , m' Jc ) 



NcgN' cg 



m Jc \D l {p c ,p' < : i 



hK'^—hR, 



\ m 'jc) ■ 



(69) 



h = Gt' (s, Q\ s') + G£(s, Q 2 , s') { [ip cv T»(p' c )] 2 -~ $ V [ip Cv r»( P ' c )f 



s P [Pc,r-(p' c ; 



F 3 = G%(s,Q 2 ,s" 



(70) 



We have taken into account that the equation G^s, Q 2 , s') = is valid in weak sense. 



5 Conclusion 

The method of construction of the electromagnetic current matrix elements for the relativistic 
two-particle composite systems with nonzero total angular momentum is developed in the 
frame of the instant form of RHD. 

The method makes use of the Wigner-Eckart theorem on the Poincare group. It enables 
one to extract from the matrix elements the reduced matrix elements - invariant form factors 
- which in the case of composite systems are generalized functions. 

The obtained current operator matrix elements satisfy the Lorentz-covariance condition 
and the conservation law. 

The modified impulse approximation (MIA) — with the physical content of the relativistic 
impulse approximation — was formulated in Refs. jU E] in terms of reduced matrix elements. 
MIA conserves Lorentz covariance of electromagnetic current and the current conservation 
law. 

The developed formalism was used in jTTj to obtain a reasonable description of the static 
moments and the electromagnetic form factors of p meson. A number of relativistic effects 
are obtained, for example, the nonzero quadrupole moment (in the case of S state) due 
to the relativistic Wigner spin rotation. In Ref . [T2*] the tenzor polarization T2o(Q 2 ) in the 
elastic e(i-scattering was calculated. The good agreement with experiments was obtained, 
especially for deuteron wave functions from Ref. JH] • These particular results argue for the 
validity of our description of the electroweak properties of composite systems with nonzero 
total angular momentum. 

This work was supported in part by the Program "Russian Universities-Basic Researches" 
(Grant No. 02.01.013) and Russian Ministry of Education (Grant E02- 3.1-34). 
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